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ABSTRACT 



To better understand the hydro-dynamic flow behavior in turbulence, Particle-Fluid flow have been studied using our 
. Direct Numerical(DNS) based software DSM on MUSCL-QUICK and finite volume algorithm. The particle flow was 
studied using Eulcrian-Eulcrian Quasi Brownian Motion(QBM) based approach. The dynamics is shown for various 
[ particle sizes which are very relevant to spray mechanism for Industrial applications and Bio medical applications. 

II. INTRODUCTION AND MOTIVATION 

V'- 

. Particle laden flows are of great interest in Chemical, Spray Industrial applications or Biomedical applications. 
rj"^ ' Knowledge of particle transport and concentration properties are crucial for experimental design of such applications. 
^ ■ Numerical simulation coupling Lagrangian tracking in discrete carrier phase of turbulence with Direct Numerical 
■ Simulation(DNS) simulation phase provides a robust tool to investigate such flows. Many applications of this type 
of flow are aerosol [Nitrogen and other post combustion particles] particle flow studies in post combustion chamber 
p^ for the design and performance optimization of aviation engine and also for spray mechanism in nano particle spray 
| Industries. In studies of arterial blood flow, these simulations help to visualize the blood flow in cardio-vascular 
Q4 region. Numerical Simulation coupling Lagrangian approach of particle tracking along with DNS (Direct Numerical 
Simulation) phase of simulation provides a powerful tool for such flow such in any geometry. 
I ■ In all Industrial applications of the boundary wall problem, boundaries are not smooth in true sense and the presence 
' of roughness cause additional energy dissipation enhancing mixing in the particle-fluid mixture. Most commercial 
C*") . software like Fluent and CFX are unable to incorporate such physics. There are two commonly used methods to 
0^ ■ simulate fluid-particle flows:Eulerian-Lagrangian and Eulerian-Eulerian, and both of these have been used to model 
the settling of particles in an in compressible fluids. The Eulerian-Lagrangian technique treats the fluid as a continuous 
medium described by Navier-Stokes equation modified to account for the fluid-particle interaction where particles are 
r*"** ' considered as point in fluid such that Newton's second law can be applied separately to each particle to track its 
. motion in a Lagrangian frame of reference. Particle-particle interaction and particle-fluid interactions are modeled 
for each particle. When particle numbers become large, particle particle interaction and turbulence modifications 
become expensive to study such kind of flows. This even sometimes becomes computationally challenging to retrieve 
data fill ]. These effects become important in the calculation which makes DNS simulation approach to become 
numerically Numerical computation based on separate Eulerian balance can provide very good alternative approach 
to such problems which are numerically not so expensive in comparison to grid size. Eulerian approach is based on 
separate balance equation for both phases through inter phase coupling terms. Such Eulerian-Eulerian DNS approach 
have been validated for the case of particles with low inertia which follow their carrier fluid almost instantaneously 
due to their small response time compared to integral time scale of turbulence [Druzini et al] [f|. In case of inertial 
particles, where particle response time scale is comparable to integral time scale, additional effects have to be taken into 
account. As pointed out by Fevrier et al Q @ calculations, particle phase transport equations must take account all 
dispersion effects due to local random motion which is induced by particle-particle interaction and particle-turbulence 
interaction. In fact, complex nonlinear particle-particle interaction exists in the vicinity of wall. 

Following Fevrier et al pj , a conditional average of all dispersed phases with respect to the carrier phase flow allows 
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the derivation of instantaneous mesoscopic particle fields and instantaneous Eulcrian balance equation which can 
be calculated by taking into account for the effect of random motion. The detail calculation of energy dissipation 
parameter also takes into account of random motion of the dispersed phase. 

Using forced isotropic turbulence simulations, Fevrier et al 0] showed that uncorrelated quasi Brownian motion of 
the particles increases with inertia(Stokes Number). 

In cases where particle relaxation time is comparable to Lagrangian integral time scale, kinetic energy of the quasi 
Brownian motion is about 30% of total kinetic energy of the dispersed phase. 

The importance of Quasi Brownian Motion(QBM) is illustrated in a preliminary test in case of decaying ho mog eneous 
isotropic turbulence. The Eulerian model is then applied to the experimental case of Snyder & Lumley [la ] which 
has been previously simulated using Lagrangian tracking approach by Elghobashi & Trusdel Q. This allows us to 
compare results of simulation with Lagrangian approach and also with the experimental data. The next section will 
illstrate the Eulerian model and related calculation. 



III. THE EULERIAN MODEL 



Eulerian equations for the dispersed phase may be derived using the methods which consists of volume filtering of 
the separate, local instantaneous equations accounting for the inter facial jump conditions [Druzini et al[f|]. Such an 
averaging approach is very restrictive because particle size and inter particle distance have to be smaller than the 
smallest turbulence scale. 

A different approach in the framework of kinetic theory is the statistical approach. In analogy to the derivation 
of Navier Stokes equation by non equilibrium statistics by Chapman & Cowling Q, a point probability density 
function(pdf) f p (l)(c p ; x pi t) defining the local instantaneous probable number of particle centers with the given 
translation velocity u p = c p , is defined. This function obeys Boltzmann type of kinetic equation which accounts for 
the momentum exchange with carrier fluid and other inter particle forces and external nonlinear force field. 
Reynold-averaged Transport equations for the first moment such as particle concentration, mean velocity and particle 
kinetic stress can be derived directly from this probability density function of Kinetic equation [Simonin et. al.[lS]]. 
To derive local instantaneous Eulerian Equation in dilute flows(without turbulence modifications of the particles), 
Fevrier et al 0] proposed an approach that uses averaging over all dispersed phases over single carrier phase condition. 
Such an averaging leads to conditional particle dispersion function(pdf) for the dispersed phase defined by, 

/ P (1) (c P ;x,t,fff) - (Wp^CpjMltff) (l) 

The W p 's are the position and velocities of corresponding particle with time. With this definition, one may define 
local instantaneous particulate velocity field defined as "Mesoscopic Eulerian Particle Velocity Field". This field 
is obtained by averaging discrete particle velocities measured at a particular position and time for all particle-flow 
realizations and at given carrier-phase realization. 

Such an averaging leads to a conditional "pdf for the dispersed phase, 

u p {u,t,H{) = — J Cp/p (1) (c p ;x,i, H f )dc p (2) 



Here 

(i) 



n 



p 



//p (1) (c p ;x,t,Ff)dcp ( 3 ) 



is the mesoscopic particle number density. Application to the conditional averaging procedure to the kinetic equation 
governing the particle pdf directly leads to the transport equation for the first moments of number density and 
mesoscopic Eulerian velocity. 



d d 



-n 



dt p dxi 



ipUp.i = (4) 



Here <5cr p "ij is the mesoscopic kinetic stress tensor of the particle Quasi Brownian velocity distribution. Our calculation 
showed that this term is non negligible for inertial particles in turbulent flow. 
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A. The Stress Tensor of Quasi Brownian Motion(QBM) 

The stress term in eqn(5) arises from the ensemble average of the non linear term in the transport equation for 
particle momentum, 

TipdoWj = / (cp -Mp,i)(cp,j - Up,j)/p (1) (c p ;x,£, H t )dc p 

J (6) 

= ™ P £up,i<5«p,i 

When Euler or Navier Stokes equation is derived from kinetic gas theory, the trace of 5u p 5u p j is interpreted as 
temperature(ignoring Boltzman constant and molecular mass) and related to pressure by equation of state. 
In Euler or Navier Stokes equation, temperature is defined as the uncorrelated part of kinetic energy. And uncorrelatcd 
part of kinetic energy is defined as, 

In analogy to Euler or Navier Stokes equation, product of uncorrelated kinetic energy and particle number density is 
defined as uncorrelated Quasi Brownian Pressure(QBP) as, 

^ P = n P ^Sq p 2 (8) 

When particle number density becomes non uniform, as in the case of compressible gas, the pressure tends to homog- 
enize particle number density. 

The non diagonal part of stress tensor can be identified in analogy to Navier Stokes Equation as viscous tcrm(0;j). 
The momentum transport equation (5) becomes, 

d d n p d ~ d 

n " P ^"p4 + "p w pj ^-"p,i = --rKi - - q^: p p + a^j.p + n p9i ( 9 ) 

Moreover, it can be shown mathematically that equation (9) without pressure like term leads to nonphysical solution. 

B. Simulation with and without QBM 

Preliminary simulations were performed without any stress term related to QBM. Particles tend to accumulate 
rapidly in smaller region causing unnecessary high particle number density. This causes numerical simulation to 
fail. In order to ensure that failure is not caused by numerical problem, different simulations with different Reynold 
Number was carried out showing same result. 

Simulations with quasi Broqnian pressure(QBP) and without quasi Brownian viscous stress were done on the above 
test cases. Fevrier et al @ measured in forced homogeneous isotropic turbulence, a mean quasi Brownian kinetic 

energy dq 2 proportional to the mean mesoscopic kinetic energy q p 2 = -(up,iW Pl i) with a proportionality coefficient 
depending on the Stokes number. 

The relation between the quasi Brownian kinetic energy and mean resolved kinetic energy is used, 

<5<7 P 2 = 5q p 2 (10) 

Such a QBP model allows all the test cases that failed without quasi Brownian stress term to simulate. Compared to 
the value obtained by Fevrier et al Q, relation (10) overestimates the quasi Brownian kinetic energy. When viscous 
stress term is incorporated in the calculation, it reduces the pressure. 

In order to quantify the effect of particle confinement, the normalized variance of particle number density is introduced, 

(n(x,t)n(x + r,t)) 
((n{x,t)) 



Figure 1 compares time evolution of g(0, t) from simulations with and without QBP. The quasi Brownian pressure is 
found to affect the particle conglomeration or segregation significantly. 
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FIG. 1: Fig 1. Normalized variance distribution of particle number density without QBP(green color) and with QBP(red color) 

Figure 2 shows kinetic energy spectra of the carrier phase and the dispersed phase with and without QBP. Our 
simulations show same behavior of simulations that without QBP case, particle kinetic energy of small scale turbulence 
becomes larger than the carrier phase in contrast to the results of Fevrier et al Q . This is probably due to nonphysically 
large particle accumulation in the carrier phase of turbulent flow which is the region of high shear strain and low 
vorticity. The vorticity development can be another useful tool to study such flows. Figure 5 shows the vorticity 
development for the carrier phase as well as dispersed phase. In presence of QBP contribution limiting particle 
segregation, particle enstrophy behavior shows similar effect. Figure 3 shows particle enstrophy distribution with and 
without QBP. The quasi Brownian viscous term plays a major role by inducing a strong small scale strong dissipation 
effect in addition to drag force. Figure 3 shows the effect by characterizing the temporal increase of enstrophy function. 
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(a)Fig 2. non dimensional energy spectra of (b)Fig 3. Enstrophy distribution in dispersed 

+ (red) carrier phase o(pink) without QBM phase with X (green) and without + QBP 

and x (blue) with QBM 



IV. PARTICLE DYNAMICS AND DISPERSION 

Particle dispersion in Lagrangian simulations is usually measured by tracking individual particle path and calcu- 
lating the variance of the relative displacement. In our simulation, Eulerian-Eulerian simulation was performed with 
one way coupling. We also studied particle dynamics in simulations with gravity and without gravity to study the 
gravitational effect on dynamics with nonlinear force. The Eulerian-Eulerian simulation was carried out in a coupled 
way so that carrier phase turbulence affects the dispersed phase through particle dynamics. The interaction in the 
momentum equation is drag force and nonlinear Basset force in the limit of large density ratio p p /pi- 
The characteristic relaxation time is computed by formulation, 



1 



<V(*)) = ^E K,i(*)-*pj(*o)r 



(12) 
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The characteristic polynomial of the Jacobian Matrix is given by, 



(«-A)(u-A)-y) 



/ 5 

(lip - A) (lip - A) 



V 
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(13) 



Particle dispersion can then be related to the time derivative of the quantity |13| 

D P L (t) = ~(X p 2 (t)) (14) 

Eulerian simulation can not track individual particle path. Particle dispersion is measured by semi-empirical method 
[l3j If the simulation is being carried out with colored particles and transport equation is written in terms of ratio of 

colored particles to total particles (c = — ). Then we can write transport equation as, 

rip 

|en p + ±an p u p = ± (15) 

Here M Pl i c is the mesoscopic velocity of colored particles. Since only the velocity of total particle is resolved, the right 
hand side term takes into account of the slip velocity between colored and mesoscopic velocity of the particle ensemble. 
Comparing to Navier-Stokes equation, this term is equivalent to molecular diffusion. Since slip velocity arises only 
from uncorrelated movement of particles, this term can be modeled as diffusion term related to quasi brownian motion. 
If the ensemble averaged mean number density fraction of colored particles {n p )C = (n p c) is uniformly stratified in 
the k direction [here in — k direction] c = C + c and fluctuations are assumed periodic with respect to the computer 
domain, then fluctuation number density of colored particles cn p can be calculated from the total colored particle 
density function. One obtains a transport equation for the fluctuation of colored particle concentration as, 

d d d d 

T^crip + T—cripUp.i = -n p Mp, k - — C + —— cri p (u p A - u p i ) (16) 
at oxi oxk oxi 

Averaging the colored number density equation (15), one obtains Reynold averaged type transport equation, 
d d d d 

dV p '° + fo Wp)C{upj)p = -— (np<5tip,i) + — (cri p (Mp,i - m p /)) (17) 
Particle dispersion term can be derived by making gradient assumption, 

«cn>p-, k ) = {npDp^^-C) (18) 
A semi empirical diffusion coefficient can be defined as, 

Dp^ = (19 ) 

This dispersion coefficient is comparable to the Lagrangian Dispersion cocfficient(13) in the long time limit of station- 
ary turbulence. Simulations without QBP likely underestimates Lagrangian dispersion. The characteristic particle 
relaxation time is computed by the formulation given by, 

Tp = imnep)n (20) 

Particle Reynold number for the drag force correction f(Re p ) is based 

f(Re p ) = 1 + 0.15i?e p °- 687do (21) 
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V. NUMERICAL SIMULATIONS AND RESULTS 



Particle dynamics and particle dispersion have been studied by experiments and Lagrangian computations. Ex- 
perimental results of Snyder & Lumley fl(| (referred as SL) is very robust test case for numerical simulation. They 
inserted particles with different incrtial properties into grid generated spatially decreasing turbulence and measured 
particle dynamics as well as particle dispersion. The test case has been computed with Lagrangian approach by 
Elghobashi & Truesdell(ET) Q . The carrier phase was taken as a temporarily decreasing homogeneous isotropic 
turbulence corresponding to the grid generated turbulence of SL. After an initial calculations for two turnover time 

(t = -7-), particles were inserted into the flow. Particle dynamics and dispersion were analyzed by ET 0] on particles 

similar to the experiments of Snyder & Lumley [16( for direct comparison with lagrangian simulation. Wc carried out 
particle simulation with Eulcrian-Eulcrian approach and comparison with the experimental results and Lagrangian 
simulation results were attempted. The simulations were performed on 100 3 grid. 

The carrier phase velocity is initialized at dimensionless time T=0 with a divergence velocity field (continuum condi- 
tion) such that the kinetic energy satisfies the spectrum Q 



E(k,0) = ^,o 3 r exp((-fc//c p ) 



(22) 



where Uf is the dimensionless rms velocity, k is the wave number and k p is wavenuber of peak energy. All wave 
number were normalized with fc m i n . In the present simulation, properties of carrier phase was validated against the 
properties of carrier phase turbulence of SL and ET. The spatial evolution of the flow in the experiment of SL is 



converted to a temporal evolution of the flow by t = -tt 



here U is the mean convection velocity in the experiment. 



Figure 4 shows dimensionless velocity square in the carrier phase in comparison with Lagrangian simulation(ET) 
and experiment (SL) results. To verify numerical simulation resolution, dissipation energy e is also compared to the 



(c)Fig 4. Squared dimensionless velocity 
Distribution-Lagrangian (x) and present( — ) 
simulation 



(d)Fig 5. Comparison of dissipation energy in 
carrier phase with Lagrangian simulation 



temporal change of kinetic energy — q 2 in Figure 5. Our result shows excellent agreement between dissipation and 

at 

kinetic energy decrease. So it can be assumed that numerical dissipation is negligible compared to viscous dissipation. 
In Figure 6, Reynold number from our simulation with Lagrangian simulation. The present simulation more rapid 
decay of turbulent Reynold number compared to the Lagrangian simulation(ET). This is because the integral length 
scale increase slowly in our simulation(Figure 7). The temporal evolution of integral time length scale shows similar 
qualitative behavior in eulerian simulation. Particles were inserted as in the Lagrangian simulation Uf' (ET) at 
nondimensional time T = 2.0 with same velocity as carrier phase when inserted into the turbulent flow. Particle prop- 
erties were analyzed in turbulence with and without gravity. When particles are inserted into gravity, they establish 
a mean terminal velocity vertically downward direction. The gravity constant was calculated as in the experiment. 
For all types of particle sizes, relative square velocity in the present simulation shows same qualitative behavior as 
in the Lagrangian simulation. Figure 8 and 9 shows the squared relative velocity distribution of the carrier phase 

Vt 

without and with the presence of gravity. The gravity constant is measured such that same ratio of — — . In all cases, 

uq 

present simulation shows same behavior as Lagrangian simulation. We have chosen temporarily decaying turbulence 
condition in our simulation as initial parameter. In Eulerian simulations, one does not have access to individual 
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bvuLlion of Carrier Phase Reynold N 



Evolution of integral length 



(e)Fig 6. Evolution of Reynold Number in carrier 
phase X (green) Lagrangian and +(red) present 
simulation 



(f)Fig 7. Evolution of integral length scale 
X (pink).. Lagrangian and +(green) ..present 
simulation 



particle paths. Particle dispersion can still be measured by a semi-empirical method |l3(. Since only the velocity 
of the total droplet number is resolved which brings another extra term on right hand side of the above equation. 
This term takes into account the slip velocity between colored particle and the mesoscopic velocity of particle ensemble. 

dispersed phase in excellent agreement. It is observed in the present simulation that when gravity is taken into 
account, particle dynamics was modified. The observed crossing trajectory effect is due to mean settling velocity 
of the particle and leads to decrease in integral time scale of fluid turbulence. Such an effect leads to an increase 
of effective particle stokes number. This also increases the relative squared velocity value with respect to the non 
settling value. Particle dispersion is also calculated for the dispersed phase. In order to compare with the carrier 
phase, equation (16) is solved for carrier phase without molecular diffusion. In the Eulerian simulation with gravity 
taking into account, particle dispersion is significantly lower than the simulation without gravity consistent with 
Csanady analysis 0). Calculation for all particles show that Eulerian simulation shows similar qualitative behavior 
as Lagrangian simulation and dispersion calculation is quantitatively lower than Lagrangian simulation. 



(g)Fig 8. Evolution relative 
vclocity(dimensionless) with gravity 



(h)Fig 9. Evolution of relative 
velocity(dimcnsionlcss) without gravity 



The Eulerian equations for the dispersed phase have been implemented on the slip velocity 

Re p = fo-W (23) 

The Eulerian mean square relative velocity differs from the Lagrangian mean square velocity by the quantity 
Su p 2 from QBM. Figure 12 shows the temporal development of the carrier phase u 2 . Balakin et al [l[ also studied 
uniformly sized sediment particle movements in viscous fluids using Eulcrian-Eulcrian simulation and compared with 
experimental results of Nicolai et. al [l4| for sedimenting particles. They corresponding governing equation is given 
by, 

%^ + VMa,) = (24) 
at 



Index m assigns the phase (liquid or solid), a is the volume fraction, p is the density and u is mean phase velocity. 
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Evolution of Dispersion Coefficient without Gravity 



Evolution of Dispersion Coefficient with Cruvny 
0.2 i i i i i i 



(i)Fig 10. Evolution of dispersion coefficient 
without gravity 



(j)Fig 11. Evolution of dispersion coefficient with 
gravity 



The momentum equation is given by, 
9(a m PmUm) 



at 



+ V(amPmUmMm) = -Q^mVP + Ct m Pm9 + V(amT m ) + M m + F A 



(25) 



In further to investigate our simulation approach, we carried out simulation for scdimcnting particle flow simulation 
similar to the experiment done by Nicolai et al. fbl ] Figure 11 shows the comparison between our simulation with 
Experimental results of Nicolai etl al. 14 and also Balakin et al |l| simulation results. The dispersed phase velocity is 



plotted against volume fraction of carrier phase to dispersed phase. Our simulation strongly agree with experimental 
data. . 
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Solid Volume Fraction 

(k)Fig 12: Time average settling velocity plotted 
as a function of solid volume fraction, red X ■ ■ ■ 
(Nicolai et al), blue X • ■ ■ balakin et al and green 
o • ■ ■ (our simulation) 



VI. CONCLUSION 



Particle dispersion is measured for the dispersed phase. In order to compare with the carrier phase, the equivalent 
of equation(16) is solved for carrier phase without molecular diffusion. A preliminary model of Quasi Brownian 
Motion was used to study unresolved kinetic particle energy. This model allowed simulations to compare with the 
experimental results of Snyder & Lumley fl6| . We also have done for sedimenting particles similar to experiment done 
by Nicolai et. al [l4| and compared with Balakin simulation results. [l[ which strongly shows that Quasi Brownian 
ensemble approach is very strong alternative tool for two phase flow modeling. 
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